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CN : ABSTRACT 
O \ We present deep near-infrared (NIR) J^, H, and i^^-band ISAAC imaging of the WFPC2 field of the 

^ \ Hubble Deep Field South (HDF-S). The 2.5' x 2.5' high Galactic latitude field was observed with the 

VLT under the best seeing conditions with integration times amounting to 33.6 hours in J a, 32.3 hours 
in H, and 35.6 hours in Kg. We reach total AB magnitudes for point sources of 26.8, 26.2, and 26.2 
respectively (3cr), which make it the deepest ground-based NIR observations to date, and the deepest 
iCg-band data in any field. The effective seeing of the coadded images is « 0'.'45 in J^, « 0'.'48 in H , and 
. « 0'.'46 in Kg. Using published WFPC2 optical data, we constructed a iiTs-limited multicolor catalog 

^ ■ containing 833 sources down to Kl"^^ < 26, of which 624 have seven-band optical-to-NIR photometry. 

, These data allow us to select normal galaxies from their rest-frame optical properties to high redshift 

■ (-2^4). The observations, data reduction and properties of the final images are discussed, and we 

address the detection and photometry procedures that were used in making the catalog. In addition, we 
present deep number counts, color distributions and photometric redshifts of the HDF-S galaxies. We 
' find that our faint i^s-band number counts are flatter than published counts in other deep fields, which 

I might reflect cosmic variations or different analysis techniques. Compared to the HDF-N, we find many 

■~-«^ . galaxies with very red V — H colors at photometric redshifts 1.95 < Zp/iot < 3.5. These galaxies are 

' bright in Kg with infrared colors redder than Jg — Kg > 2.3 (in Johnson magnitudes). Because they are 

extremely faint in the observed optical, they would be missed by ultraviolet-optical selection techniques, 
Q ■ such as the U-dropout method. 

Subject headings: cosmology: observation — galaxies: evolution — galaxies: high-redshift — galaxies: 
^ ' photometry 
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1. INTRODUCTION history, stellar population ages, and masses are not well 

^ In the past decade, our ability to routinely identify and ^'^f'!)^'^' • ; •, • , -r,, ,; . , ; ;• , 

, 11 i J J- i i 1-1. J 4-- 11 J More importantly, it is unclear if the ultraviolet-optical 

^ 1 systematically study distant galaxies has dramatically ad- , ,. , , . , •,, . r • r 

^ , J 1 1 J r , • 1 1 1 -r; • T selection technique alone will give us a fair census of the 

■ - ' vanced our knowledge oi the high-redshift universe, in , , ; • , o • i ■ . i 

,-1 a; • i TT J 4. 4. I, • -J 1 i galaxy population at z ^ 3 as it requires galaxies to have 

particular, the ethcient U-dropout technique (Steidel et f - i n i. ■ i . r i • i ; ^ . 

al. 1996a,b) has enabled the selection of distant galaxies ^^S^. far-ultraviolet surface brightnesses due to on-gomg 

r 4-- 1 ■ ■ ■ ■ 1 u 4- i ■ ■ spatially compact and relatively unobscured massive star 

from optical imaging surveys using simple photometric cri- / , .■' ^ , i , , i • , i • i i i ^ 

. . ,1 innn 4- 4-^, t u i i • formation. We know that there exist highly obscured 

teria. Now more than 1000 of these Lyman break galaxies , . , ; , , . , , ,. -i ^ 

/T lar^ \ , -11 c ji-^o ju galaxies, detected m sub-mm and radio surveys (Smail et 

(LBGs) are spectroscopically confarmed at z > 2, and have °, „r.rsr,\ i . • n n • ; i i /V, • , 

^ \_- . . . ,1 ,1- ,-1 ?^ , • //-I- al. 2000), and optically famt hard A-ray sources Gowie et 

been subiect to targeted studies on spatial clustering (Gi- , „„„/V> . / nr.r.^\ , i - i i i -n ;i ; ii 

T p T-,- 1 • -4- 11- 4-- /T) 4-4-- • al- 2001; Barger et al. 2001 at high redshift that would 

avalisco & Dickinson 2001), internal kinematics (Pcttini , ' °, t i i • i i 

et al. 1998, 2001), dust properties (Adelberger & Steidel ^ ^^^^^^f^ ^t'J;? .'^ u ' densities are 

2000), and stellar composition (Shapley et al. 2001; Pa- low compared to LBGs and they might represent rare pop- 

. 1 -p.. 1 . p ^ r,nn1^ A 14.1, UT-D/"! ulations or transient evolutionary phases. In addition, the 

povich, Dickinson, & l<erguson 2001). Although LBGs are . ^ , , „. ;• , , . , , ■ , 

4.1 1 4_ , J. i 1 r J- 4. 4. 1 • 4. J 4. majority of present-day elliptical and spiral galaxies, when 

among the best studied classes of distant galaxies to date, i i r, , , ; • n r ;? 

r . 4-- 1-1 4-1, • • 4- f 4-- placed at z ~ 3, would not satisfy any of the current se- 

many of their properties like their prior star formation ' j j 

^ Based on service mode observations collected at the European Southern Observatory, Paranal, Chile (ESO Programme 164.O-0612). Also 
based on observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope Science Institute which is operated by 
AURA, Inc., under NASA contract NAS5-26555. 

2 Leiden Observatory, P.O. Box 9513, NL-2300 RA, Leiden, The Netherlands 
^ Max-Plank-Institut fiir Astrophysik, P.O. Box 1317, D-85741, Garching, Germany 
^ Max-Plank-Institut fiir Astronomic, D-69117, Heidelberg, Germany 
^ European Southern Observatory, D-85748, Garching, Germany 
^ California Institute of Technology, MS 105-24, Pasadena CA 91125, USA 
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lection techniques for high-redshift galaxies. Specifically, 
they would not be selected as U-dropout galaxies because 
they are too faint in the rest-frame UV. It is much easier 
to detect such galaxies in the near- infrared (NIR), where 
one can access their rest-frame optical light. Furthermore, 
observations in the near-infrared allow the comparison of 
galaxies of different epochs at fixed rest-frame wavelengths 
where long-lived stars may dominate the integrated light. 
Compared to the rest-frame far-UV, the rest-frame opti- 
cal light is less sensitive! to the effects of dust extinction 
and on-going star formation, and provides a better tracer 
of stellar mass. By selecting galaxies in the near-infrared 
Ks-hand, we expect to obtain a more complete census of 
the galaxies that dominate the stellar mass density in the 
high-redshift universe, thus tracing the build-up of stellar 
mass directly. 

In this context we initiated the Faint InfraRed Ex- 

tragalactic Survey (FIRES; Franx ct al. 2000), a large 
public program carried out at the Very Large Telescope 
(VLT) consisting of very deep NIR imaging of two se- 
lected fields. We observed fields with existing deep optical 
WFPC2 imaging from the Hubble Space Telescope (HST): 
the WPFC2-ficld of Hubble Deep Field South (HDF-S), 
and the field around the z « 0.83 cluster MS1054-03. The 
addition of NIR data to the optical photometry is required 
not only to access the rest-frame optical, but also to deter- 
mine the redshifts of faint galaxies from their broadband 
photometry alone. While it may be possible to go to even 
redder wavelengths from the ground, the gain in terms of 
effective wavelength leverage is less dramatic compared to 
the threefold increase going from the I to fiT-band. This is 
because the Xg-band is currently the reddest band where 
achievable sensitivity and resolution are reasonably com- 
parable to deep space-based optical data. Preliminary re- 
sults from this program were presented by Rudnick et al. 
(2001, hereafter ROl). 

Here we present the full NIR data set of the HDF-S, to- 
gether with a iiTs-selected multicolor catalog of sources in 
the HDF-S with seven-band optical-to-infrared photome- 
try (covering 0.3 — 2.2/xm), unique in its image quality and 
depth. This paper focusses on the observations, data re- 
duction and characteristic properties of the final images. 
We also describe the source detection and photometric 
measurement procedures and lay out the contents of the 
catalog, concluding with the NIR number counts, color 
distributions of sources, and their photometric redshifts. 
The results of the MS1054-03 field will be presented by 
Forster Schreiber et al. (2002) and a more detailed expla- 
nation of the photometric redshift technique can be found 
in Rudnick et al. (2001, 2002b). Throughout this paper, all 
magnitudes are expressed in the AB photometric system 
(Oke 1971) unless explicitly stated otherwise. 

2. OBSERVATIONS 

2.1. Field Selection and Observing Strategy 

The high Galactic latitude field of the HDF-S is a natu- 
ral choice for follow-up in the near-infrared given the exist- 
ing ultradeep WFPC2 data in four optical filters (Williams 
et al. 1996, 2000; Casertano et al. 2000). The Hubble 
Deep Fields (North and South) are specifically aimed at 
constraining cosmology and galaxy evolution models, and 
in these studies it is crucial to access rest-frame optical 



wavelengths at high redshift through deep infrared obser- 
vations. Available ground-based NIR data from SOFI on 
the NTT (da Costa et al. 1998) are not deep enough to 
match the space-based data. To fully take advantage of 
the deep optical data requires extremely deep wide-field 
imaging in the infrared at the best possible image qual- 
ity; a combination that in the southern hemisphere can 
only be delivered by the Infrared Spectrometer And Ar- 
ray Camera (ISAAC; Moorwood 1997), mounted on the 
Nasmyth-B focus of the 8.2 meter VLT Antu telescope. 
The infrared camera has a 2.5' x 2.5' field of view similar 
to that of the WFPC2 (2.7' x 2.7'). ISAAC is equipped 
with a Rockwell Hawaii 1024 x 1024 HgCdTe array, offer- 
ing imaging with a pixel scale of 0'.'147 pix~^ in various 
broad and narrow band filters. 

Our NIR imaging consists of a single ISAAC point- 
ing centered on the WFPC2 main-field of the HDF-S 
(a = 22''32'"55.464, S = -60°33'05.01", J2000) in the J^, 
H and Ks filters, which gives good sampling of rest-frame 
optical wavelengths over the redshift range 1 < z < 4. 
The Jg filter is being established as the new standard 
broadband filter at w 1.24/xm by most major observato- 
ries (Keck, Gemini, Subaru, ESO), and is photometrically 
more accurate than the classical J because it is not cut off 
by atmospheric absorption. It is a top-hat filter with sharp 
edges, practically the same effective wavelength as the nor- 
mal J filter, and half-transmittance points at 1.17/um and 
1.33/J,m. We used the filter which is bluer and nar- 
rower than standard K, but gives a better signal-to-noise 
ratio (SNR) for faint sources because it is less affected 
by the high thermal background of the atmosphere and 
the telescope. The ISAAC H and Kg filters are close to 
those used to establish the faint IR standard star system 
(Persson et al. 1998), while the Jg filter requires a small 
color correction. The WFPC2 filters that are used are 
F300W, F450T4^, i^606T4^ and F814T4^ which we will call 
U^ao, -B450, Veo6 and /814, respectively, where the subscript 
indicates the central wavelength in nanometers. 

The observing strategy for the HDF-S follows estab- 
lished procedures for ground-based NIR imaging. The 
dominance of the sky background and its rapid variability 
in the infrared requires dithering of many short exposures. 
We used a 20" jitter box in which the telescope is moved in 
a random pattern of Poissonian offsets between successive 
exposures. This jitter size is a trade-off between keeping 
a large area at maximum depth and ensuring that each 
pixel has sufficient exposures on sky. Individual exposures 
have integration times of 6 x 30 s in Jg, 6 x 20 s in i?, 
and 6 x 10 s in fCs (subintegrations x detector integration 
times). We requested service mode observations amount- 
ing to 32 hours in each band with a seeing requirement 
of < 0'.'5, seeing conditions that are only available 25% 
of the time at Paranal. The observations were grouped 
in 112 observation blocks (OBs), each of which uniquely 
defines a single observation of a target, including point- 
ing, number of exposures in a sequence, and filter. The 
calibration plan for ISAAC provides the necessary cali- 
bration measurements for such blocks, including twilight 
flats, detector darks, and nightly zero points by observ- 
ing LCO/Palomar NICMOS standard stars (Persson et al. 
1998). 
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Fig. 1. — Shown arc the raw data in the filters Js {dotted line or circles), H (dashed line or squares), and Kg {solid line or triangles), 
(a) Histogram of the median seeing in tlie raw ISAAC images weighted by the weight function of Eq. 2 used to combine the images, (b) 
Relative instrumental counts in a 3" radius aperture of four bright non-saturated stars in individual sky-subtracted exposures, plotted 
against Julian Date. The relative increase in counts, slightly dependent on wavelength, after cleaning and re-aluminization of the mirror 
directly reflects the increase in efHciency of the telescope, because the sky background levels (c) remained the same. Presumably, the photons 
were scattered by dirt rather than absorbed before cleaning, (d) Nightly sky variations are largest and most rapid in the i?-band and mean 
sky levels are highest at the beginning and ending of the night. Js-band varies less and peaks at the start of the night, whereas Jf-band levels 
are most stable. 



2.2. Observations 

The HDF-S was observed from October to December 
1999 and from April to October 2000 under ESO program 
identification 164.O-0612(A). A summary of the observa- 
tions is shown in 1. We obtained a total of 33.6, 32.3 and 
35.6 hours in Js, H and Kg, distributed over 33, 34 and 55 
OBs, or 1007, 968 and 2136 frames, respectively. This rep- 
resents all usable data, including aborted and re-executed 
OBs that were outside weather specifications or seeing con- 
straints. In the reduction process these data arc included 
with appropriate weighting (see section 3.4). Sixty-eight 
percent of the data was obtained under photometric con- 
ditions and the average airmass of all data was 1.25. A 
detailed summary of observational parameters with point- 
ing, observation date, image quality and photometric con- 
ditions can be found on the FIRES homepage on the World 
Wide Web (http://www.strw. leidenuniv-nl/'fires). 

An analysis of various observational parameters reveals 



some surprising trends in the data, whereas other expected 
relations are less apparent. An overview is given in Figure 
1. The median seeing on the raw images is better than 0'.'5 
in all bands, with the seeing of 90% of the images in the 
range 0'.'4 — 0'.'65, as can be seen in Figure 1. 

Seeing may vary strongly on short timescales but it is 
not related to any other parameter. The most drastic 
trend in the raw data is the change of sensitivity with 
date. Since the cleaning and re-ahimization of the pri- 
mary mirror in March 2000 the count rates of bright stars 
within a « 3" aperture increased by -|-29% in Jg, -|-45% in 
H, and +45% Kg, which is reflected by a change in zero 
points before and after this date. Because the average NIR 
sky levels remained the same, this increase proportionally 
improved the achievable signal-to-noise for background- 
limited sources. The change in throughput was caused by 
light scattering, which explains why the sky level remained 
constant. Sky levels in Jg and H, dominated by airglow 
from OH-emission lines in the upper atmosphere (typi- 
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cally 90 km altitude) , vary unpredictably on the timescale 
of minutes, but also systematically with observed hour. 
The average sky level is highest at the beginning and end 
of each night with peak-to-peak amplitudes of the vari- 
ation being 50% relative to the average sky brightnesses 
over the night. The background in is dominated by 
thermal emission of the telescope, instrument, and atmo- 
sphere and is mainly a function of temperature. The Ks 
background is the most stable of all NIR bands and only 
weakly correlated with airmass: our data do not show a 
strong thermal atmospheric contribution, which should be 
proportional to atmospheric path length. We take into ac- 
count the variations of the background and seeing through 
weighting in the data coadding process. 

3. DATA REDUCTION 

The reduction process included the following steps: 
quality verification, flat-fielding, bad pixel correction, sky 
subtraction, distortion correction, registration, photomet- 
ric calibration and weighting of individual frames, and 
combination into a single frame. We used a modified 
version of the DIMSUM'' package and standard routines 
in IRAF^ for sky subtraction and coadding, and the 
ECLIPSE^ package for creating the flatfields and the ini- 
tial bad pixel masks. We reduced the ISAAC observations 
several times with an increasing level of sophistication, ap- 
plying corrections to remove instrumental features, scat- 
tered light, or clear artifacts when required. Here we de- 
scribe the first version of the reduction (vl.G) and the last 
version (v3.0), leaving out the intermediate trial versions. 
The last version produced the final Jg, H and Kg images, 
on which the photometry (see section 5) and analysis (see 
section 8) is based. 

3.1. Flatfields and Photometric Calibration 

We constructed flatfields from images of the sky taken 
at dusk or dawn, grouped per night and in the relevant 
filters, using the flat routine in ECLIPSE, which also 
provided the bad pixel maps. We excluded a few flats of 
poor quality and flats that exhibited a large jump between 
the top row of the lower and bottom row of the upper half 
of the array, possibly caused by the varying bias levels of 
the Hawaii detector. We averaged the remaining nightly 
flats per month, and applied these to the individual frames 
of the OBs taken in the same month. If no flatfield was 
available for a given month we used an average flat of 
all months. The stability of these monthly flats is very 
good and the structure changes little and in a gradual 
way. We estimate the relative accuracy to be 0.2 — 0.4% 
per pixel from the pixel-to-pixel rms variation between dif- 
ferent monthly flats. Large scale gradients in the monthly 
flats do not exceed 2%. We checked that standard stars, 
which were observed at various locations on the detector, 
were consistent within the error after flatfielding. 

Standard stars in the LCO/Palomar NICMOS Ust (Pers- 
son et al. 1998) were observed each night, in a wide five- 
point jitter pattern. For each star, on each night, and in 

DIMSUM is the Deep Infrared Mosaicing Software package developed by Peter Eisenhardt, Mark Dickinson, Adam Stanford, and John Ward, 
and is available via ftp to ftp://iraf.noao.edu/contrib/dimsumV2/ 

^ IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the AURA, Inc., under cooperative agreement 
with the NSF. 

^ ECLIPSE is an image processing package written by N. Devillard, and is available at ftp://ftp.hq.eso.org/pub/eclipse/ 



each filter, we measured the instrumental counts in a cir- 
cular aperture of radius 20 pixel (2'.'94) and derived zero 
points per night from the magnitude of that star in the 
NICMOS list. We identify non-photometric nights after 
comparison with the median of the zero points over all 
nights before and after re-aluminization in March 2000 
(see section 2.2). The photometric zero points exhibit a 
large increase after March 2000 but, apart from this, the 
night-to-night scatter is approximately 2%. We adopted 
the ni(;aii of the zero points after March 2000 as our ref- 
erence value. See Table 2 for the list of the adopted zero 
points. By applying the nightly zero points to 4 bright 
unsaturated stars in the HDF-S, observed on the same 
night under photometric conditions, we obtain calibrated 
stellar magnitudes with a night-to-night rms variation of 
only « 1 — 1.5%. No corrections for atmospheric absorp- 
tion were required because the majority of the science data 
wcirc obtained at similar airmass as the standard star ob- 
servations. In addition, instrumental count rates of HDF-S 
stars in individual observation blocks reveal no correlation 
with airmass. We used the calibrated magnitudes of the 
4 reference stars, averaged over all photometric nights, to 
calibrate every individual exposure of the photometric and 
non-photometric OBs. The detector non-linearity, as de- 
scribed by Amico et al. (2001), affects the photometric 
calibration by < 1% in the if-band, where the exposure 
levels are highest. Because the effect is so small, we do not 
correct for this. We did not account for color terms due 
to differences between the ISAAC and standard filter sys- 
tems. Amico et al. (2001) report that the ISAAC H and 
Kg filters match very well those used to establish the faint 
IR standard star system of Persson et al. (1998). Only the 
ISAAC Js filter is slightly redder than Persson's J and this 
may introduce a small color term, « —0.04 • (J — K)lco- 
However, the theoretical transformation between ISAAC 
magnitudes and those of LCO/Palomar have never been 
experimentally verified. Furthermore, the predicted color 
correction is small and could not be reproduced with our 
data. In the absence of a better calibration we chose not to 
apply any color correction. We did apply galactic extinc- 
tion correction when deriving the photometric redshifts, 
see section 6, but it is not applied to the catalog. 

As a photometric sanity check, we compared 2" circu- 
lar diameter aperture magnitudes of the brightest stars in 
the final (version 3.0, described below) images to magni- 
tudes based on a small fraction of the data presented by 
ROl. Each data set was independently reduced, the cal- 
ibration based on different standard stars, and the shal- 
lower data were obtained before re-aluminization of the 
primary mirror. The magnitudes of the brightest sources 
in all bands agree within 1% between the versions, indi- 
cating that the internal photometric systematics are well 
under control. For the NIR data, the adopted transfor- 
mations from the Johnson (1966) system to the AB sys- 
tem are taken from Bcssell & Brett (1988) and we ap- 

Pfy Js,AB = Js,Vega + 0.90, Hab = Hyega + 1-38 and 
Ks,AB = Ks,Vega + 1-86. 
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3.2. Sky Subtraction and Cosmic Ray Removal 

The rapidly varying sky, typically 25 thousand times 
brighter than the sources we aim to detect, is the pri- 
mary limiting factor in deep NIR imaging. In the longest 
integrations, small errors in sky subtraction can severely 
diminish the achievable depth and affect faint source pho- 
tometry. The IRAF package DIMSUM provides a two- 
pass routine to optimally separate sky and astronomi- 
cal signal in the dithered images. We modified it to 
enable handling of large amounts of data and replaced 
its co-adding subroutine, which assumes that the im- 
ages are undersampled, by the standard IRAF task IM- 
AGES.IMMATCH.IMCOMBINE. The following is a brief 
summary of the steps performed by the REDUCE task in 
DIMSUM. 

For every science image in a given OB a sky image is 
constructed. After scaling the exposures to a common me- 
dian, the sky is determined at each pixel position from a 
maximum of 8 and a minimum of 3 adjacent frames in 
time. The lowest and highest values are rejected and the 
average of the remainder is taken as the sky value. These 
values are subtracted from the scaled image to create a 
sky subtracted image. A set of stars is then used to com- 
pute relative shifts, and the images are integer registered 
and averaged to produce an intermediate image. All astro- 
nomical sources are identified and a corresponding object 
mask is created. This mask is used in a second pass of sky 
subtraction where pixels covered by objects are excluded 
from the estimate of the sky. 

The images show low-level pattern due to bias varia- 
tions. Because they generally reproduce they are auto- 
matically removed in the skysubtraction step. We find 
cosmic rays with DIMSUM, using a simple threshold al- 
gorithm and replacing them by the local median, unless a 
pixel is found to have cosmic rays in more than frames per 
OB. In this case the pixel is added to the bad pixel map 
for that OB. 

3.3. First Version and Quality Verificaiion 

The goal of the first reduction of the data set is to pro- 
vide a non-optimized image, which we use to validate and 
to assess the improvements from more sophisticated im- 
age processing. The first version consists of registration 
on integer pixels and combination of the sky subtracted 
exposures per OB. For each of the 122 OBs, we created an 
average and a median combined image to verify that cos- 
mic rays and other outliers were removed correctly, and 
we visually inspected all 4149 individual sky subtracted 
frames as well, finding that many required further pro- 
cessing as described in the following section. Finally, we 
generated the version 1.0 images (the first reduction of the 
full data set) by integer pixel shifting all OBs to a com- 
mon reference frame, and coaveraging them into the Jg, H 
and Ks images. While this first reduction is not optimal 
in terms of depth and image quality, it is robust owing to 
its straightforward reduction procedure. 

3.4. Additional Processing and Improvements 

The individual sky subtracted frames are affected by a 
number of problems or instrumental features, which we 

ISAAC home page: http;//www.eso.org/instruments/isacLC 

We have noticed that the mapping solution changed slightly after 



briefiy describe below, together with the applied solutions 
and additional improvements that lead to the version 3.0 
images. The most important problems are: 

• Detector bias residuals, most pronounced at the 
rows where the read-out of the detector starts 
at the bottom (rows 1,2,...) and halfway (rows 
513, 514, ...), caused by the complex bias behaviour 
of the Rockwell Hawaii array. These variations are 
uniform along rows, and we removed the residual 
bias by subtracting the median along rows in 
individual sky subtracted exposures, after masking 
all sources. 

• Imperfect sky subtraction, caused by stray light 

or rapid background variations. Strong variations 
in the backgrounds, reflection from high cirrus, 
reflected moonlight in the ISAAC optics or 
patterns of less obvious origin can lead to large 
scale residuals in the sky subtraction, particularly 
in Jg and H. For some OBs, we succesfuUy 
removed the residual patterns by splitting the 
sequence in two (in case of a sudden appearance of 
stray light), or subtracting a two-piece cubic spline 
fit along rows and columns to the background in 
individual frames, after masking all sources. We 
rejected a few frames, or masked the affected areas, 
if this simple solution did not work. 

• Unidentified cosmic rays or bad pixels. A small 
amount of bad pixels were not detected by 
ECLIPSE or DIMSUM routines but need to be 
identified because we average the final images 
without additional clipping or rejection. By 
combining the sky subtracted frames in a given 
OB without shifts and with the sources masked, 
we identified remaining cosmic rays or outliers 
through sigma clipping. We added 60 — 100 
pixels per OB to the corresponding bad pixel map. 

Several steps were taken to improve the quality and lim- 
iting depths of the version 1.0 images, the most important 
of which are: 

• Distortion correction of the individual frames and 
direct registration to the 3x3 blocked 78i4 image 

(0.119"pixcl^^), our preferred frame of reference. 
We obtained the geometrical distortion coefficients 
for the 3rd order polynomial solution from the 
ISAAC WWW-pagei°. The transformation 
procedure involves distortion correcting the 
ISAAC images, adjusting the frame-to-frame 
shifts, and finding the linear transformation to 
the WFPC2 7814 frame of reference. This linear 
transformation is the best fit mapping of source 
positions in the blocked WFPC2 /814 image to the 
corresponding positions in the corrected J^-band 
image^^. Compared to version 1.0 described in 
the previous section, this procedure increases 
registration accuracy and image quality, decreases 

remount of ISAAC in March 2000, implying a 0.1% scale difference. 
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image smearing at the edges introduced by the 
jittering and differential distortion. Given the 
small amplitude of the ISAAC field distortions, the 
eflFect on photometry is negligible. In the linear 
transformation and distortion correction step 
the image is resampled once using a third-order 
polynomial interpolation, with a minimal effect of 
the interpolant on the noise properties. 

• Weighting of the images. We substantially 
improved the final image depth and quality by 
assigning weights to individual frames that take 
into account changes in seeing, sky transparency, 
and background noise. Two schemes were applied: 
one that optimizes the signal-to-noise ratio (SNR) 
within an aperture of the size of the seeing disk, 
and one that optimizes the SNR per pixel. The 
first improves the detection efficiency of point 
sources, the other optimizes the surface brightness 
photometry. The weights Wi of the frames arc 
proportional to either the inverse scaled variance 
zpscalei X vari within a seeing disk of size Si, or 
to the inverse scaled variance per pixel, where the 
scaling zpscalci is the flux calibration applied to 
bring the instrumental counts of our four reference 
stars in the HDF-S to the calibrated magnitude. 

Wi,point oc {zpscalei x van x (1) 

m,extended OC (zpscale^ X worj)"^ (2) 

3.5. Final Version and Post Processing 

The final combined Jg, H, and Kg images (version 3.0) 
were constructed from the individually registered, distor- 
tion corrected, weighted and uncUpped average of the 
1007, 968, and 2136 NIR frames respectively. Ultimately, 
less than 3% of individual frames were excluded in the fi- 
nal images because of poor quality. In this step we also 
generated the weight maps, which contain the weighted ex- 
posure time per pixel. We produced three versions of the 
images, one with optimized weights for point sources, one 
with optimized weights for surface brightness, and one con- 
sisting of the best quartile seeing fraction of all exposures, 
also optimized for point sources. The weighting has im- 
proved the image quality by 10 15% and the background 
noise by 5-10%, and distortion correction resulted in sub- 
pixel registration accuracy between the NIR images and 
/814-band image over the entire field of view. 

The sky subtraction routine in DIMSUM and our addi- 
tional fitting of rows and columns (sec section 3.4) have in- 
troduced small negative biases in combined images, caused 
by systematic oversubtraction of the sky which was skewed 
by light of the faint extended PSF wings or very faint 
sources, undetectable in a single OB. Because of this, the 
flatness of the sky on large scales was limited to about 
10~^. The negative bias was visible as clearly defined or- 
thogonal stripes at P.A.w 6°, as well as dark areas around 
the brightest stars or in the crowded parts of the images. 
To solve this, we rotated a copy of the final images back 
to the orientation in which we performed sky substraction, 
fitted a 3-piece cubic spline to the background along rows 
and columns (masking all sources), re-rotated the fit, and 



subtracted it. The sky in the final images is flat to a few 
X 10"*^ on large (> 20") scales. 

4. FINAL IMAGES 

The reduced NIR J,, H, and K,^ imagers and weight 
maps can be obtained from the FIRES- WWW homepage 
(http://www.strw.leidenuniv.nl/~fires). Throughout the 
rest of the paper we will only consider the images op- 
timized for point source detection which we will use to 
assemble the catalog of sources. 

4.1. Properties 

The pixel size in the NIR images equals that of the 3x3 
blocked WFPC2 /gu-band image at 0.119" pixel'^ The 
combined ISAAC images are aligned with the HST version 
2 images (Casertano et al. 2000) with North up, and are 
normalized to instrumental counts per second. The images 
are shallower near the edges of the covered area because 
they received less exposure time in the dithering process, 
which is reflected in the weight map containing the fraction 
of total exposure time per pixel. The area of the ISAAC 
i^s-band image with weight per pixel wk > 0.95, 0.2, and 
0.01 covers 4.5,7.2 and 8.3 arcmin^, while the area used 
for our preferred quality cut for photometry (w > 0.2 in all 
seven bands) is 4.7 arcmin^. The NIR images have been 
trimmed where the relative exposure time per pixel is less 
than 1%. 

Figure 2 shows the noise-equalized _ft',,-detcction image 
obtained by division with the square root of the exposure- 
time map. The richness in faint compact sources and the 
flatness of the background are readily visible. Figure 3 
shows a RGB color composite image of the /814, Jg, and 
Ks images. The PSF of the space-based /814 image has 
been matched to that of the NIR images at FWHM« 0'.'46 
(see section 4.2) and three adjacent WFPC2 Isi4 flanking 
fields have been included for visual purposes. We set the 
linear stretch of both images to favor faint objects. Im- 
mediately striking is the rich variety in optical-NIR colors, 
even for the faint objects, indicating that the NIR obser- 
vations are very deep and that there is a wide range of 
observed spectral shapes, which can result from different 
types of galaxies over a broad redshift range. 

4.2. Image Quality 

The NIR PSF is stable and symmetric over the field with 

a gaussian core profile and an average ellipticity < 0.05 
over the Js, H, and Kg images. The median FWHM of 
the profiles of ten selected isolated bright stars is 0"45 in 
Js, 0'.'48 H, and 0'.'46 with 0'.'04 amphtude variation 
over the images. 

For consistent photometry in all bands we convolved 
the measurement images to a common PSF, correspond- 
ing to that of the 7J-band which had worst effective seeing 
(FWHM = 0'.'48). The similarity of PSF structure across 
the NIR images allowed simple gaussian smoothing for a 
near perfect match. The complex PSF structure of the 
WFPC2 requires convolving with a special kernel, which 
we constructed by deconvolving an average image of bright 
isolated non-saturated stars in the ff-band with the aver- 
age /814-band image of the same stars. Division of the 
stellar growth curves of the convolved images by the H- 
band growth curve shows that the fractional enclosed flux 
agrees to within 3% at radii r > 0'.'35. 
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Fig. 2. — The HDF-S field in the ISAAC Ks-hand divided by the square root of the weight map (based on the fractional exposure time per 
pixel) and displayed at linear scaling. The total integration time is 35.6 hours, the stellar FWHMRi 0'.'46 and the total field size is 2.85' x 2.85'. 



4.3. Astrometry 

The relative registration between ISAAC and WFPC2 
images needs to be very precise, preferably a fraction of 
an original ISAAC pixel over the whole field of view, to 
allow correct cross-identification of sources, accurate color 
information and morphological comparison between differ- 
ent bands. To verify our mapping of ISAAC to WFPC2 
coordinates, we measured the positions of the 20 brightest 
stars and compact somces in all registered ISAAC expo- 
sures, and we compared their positions with those in the 
/814 image. The rms variation in position of individual 



sources is about 0.2 — 0.3 pixel at 0.119" pixel"^ (25 — 35 
mas), but for some sources systematic offsets between the 
NIR and the optical up to 0.85 pixel (100 mas) remain. 
The origin of the residuals is unclear and we cannot fit 
them with low order polynomials. They could be real, in- 
trinsic to the sources, or due to systematic errors in the 
field distortion correction of ISAAC or WFPC2^'^. How- 
ever, for all our purposes, the effect of positional errors of 
this amplitude is unimportant. The error in absolute as- 
trometry of the HST HDF-S coordinate system, estimated 
to be less than 40 mas, is dominated by the systematic un- 
certainty in the positions of four reference stars (Casertano 



The ISAAC field distortion might have changed over the years, but this cannot be checked because recent distortion measurements are 
unavailable. The worst case errors of relative positions ax;ross the four WFPC2 chips can be O^'l (Voit et al. 1997), but is expected to be smaller 
for the HDF-S images. 
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Fig. 3. — Three-color composite image of the ISAAC field on top of the WFPC2 main-field and parts of three WFPC2 flanking fields. 
The main-field is outlined in white and North is up. The images are registered and smoothed to a common seeing of FWHMRi 0"46, coding 
WFPC2 /814 in blue, ISAAC Js in green and ISAAC Ks in red. There is a striking variety in optical-to-infrared colors, especially for fainter 
objects. A number of sources with red colors have photometric redshifts z > 2 and they are candidates for relatively massive, evolved galaxies. 
These galaxies would not be selected by the U-dropout technique because they are too faint in the observer's optical. 



et al. 2000; Williams et al. 2000). 

4.4. Backgrounds and Limiting Depths 

The noise properties of the raw individual ISAAC im- 
ages are well described by the variance of the signal col- 
lected in each pixel since both Poisson and read noise are 
uncorrelated. However, image processing, registration and 
combination have introduced correlations between neigh- 
bouring pixels and small errors in the background subtrac- 
tion may also contribute to the noise. Understanding the 
noise properties well is crucial because limiting depths and 
photometric uncertainties rely on them. 

Instead of a formal description based on the analysis of 



the covariance of correlated pixel pairs, we followed an em- 
pirical approach where we fit the dependence of the rms 
background variation in the image as a function of linear 
size N = y/A of apertures with area A. Directly measuring 
the effective flux variations in apertures of different sizes 
provides a more realistic estimate of signal variations than 
formal Gaussian scaling (j{N) = Na of the pixel-to-pixel 
noise a, as is often done. 

We measured fluxes in 1200 non-overlapping circular 
apertures randomly placed on the registered convolved im- 
ages, which were also used for photometry. We excluded all 
pixels belonging to sources detectable in Ks at the 5a level 
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Fig. 4. — Scaling relation of the measured background rins noise as a function of linear size N = VA of apertures with area A. (a) Gaussians 
are fitted to histograms of Ks counts in randomly placed apertures of increasing size, excluding pixels belonging to sources. This correctly 
accounts for pixel-to-pixel correlations and other effects, allowing us to measure the true rms variation as a function of linear size of aperture, 
(b) The A's-band results {solid points), together with the best-fit scaling relation of Eq. 3 {solid line), show that the measured variation in 
large apertures exceeds the variation expected from linear (Gaussian) scaling of the pixel-to-pixel noise {dashed line), likely due to large scale 
correlated fluctuations of the background. 



(sec section 5.1 for detection criteria). We used identical 
aperture positions for each band i and measured fluxes for 
circular aperture diameters ranging from 0'.'5 to 3". Then 
we obtained the flux dispersions by fitting a Gaussian dis- 
tribution to the histogram of fluxes at each aperture size. 
Finally, we fitted a parameterized function of linear size to 
the different dispersions: 

a,{N) = Na,{a, + b,N)/y/Wi (3) 

This equation describes the signal variation versus aper- 
ture size N over the entire image, taking into account spa- 
tial variations as a result of relative weight Wi for each 
passband i. As can be seen in Figure 4, it provides a good 
fit to the noise characteristics. The noise is signiflcantly 
higher than expected from uncorrelated (Gaussian) noise, 
indicated by a dashed line in Figure 4b. Table 3 shows the 
best fit values in all bands and the corresponding limiting 
depths. The parameter a reflects the correlations of neigh- 
bouring pixels (a > 1), which is important in the WFPC2 



images because of heavy smoothing, but also in the ISAAC 
images given the resampling from 0'.'147 to 0'.'119 pixel^^. 
The parameter b accounts for large scale correlated varia- 
tions in the background (6 > 0). This may be caused by 
the presence of sources at very faint flux levels (confusion 
noise) or instrumental features. Typically, the large scale 
correlated contribution per pixel is only 3-15% relative to 
the gaussian rms variation, but due to the N"^ proportion- 
ality the contribution to the variation in large apc^rtures 
increases to significant levels. While the signal variations 
grow faster with area than expected from a Gaussian, at 
any specific scale the variation is consistent with a pure 
Gaussian. 

From the analysis of the scaling relation of simulated 
colors we find that part of the large scale irregularities in 
the background are spatially correlated between bands. In 
particular, we measured the rms variation of the /8i4 — Vgoe 
colors directly by subtracting in registered apertures the 
/814-band fluxes from the Veoe fluxes and fltting the dis- 
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persion of the difference at each Unear size. On large scales 
rms variations are 30% smaller than predicted from Eq. 3 
if the noise were uncorrelated. Yet, if wc subtract the two 
fluxes in random apertures, the scaling of the background 
variation is consistent with the prediction. A similar effect 
is seen for the /8i4 — Js color, but at a smaller amplitude. 
The spatial coherence of the background variations be- 
tween filters and across cameras suggests that part of the 
background fluctuations may be associated with sources at 
very faint flux levels. Other contributions arc likely simi- 
lar flatflelding or skysubtraction residuals from one band 
to another. 

5. SOURCE DETECTION AND PHOTOMETRY 

The detection of sources at very faint magnitudes 
against a noisy background forces us to trade off complete- 
ness and reliability. A very low detection threshold may 
generate the most complete catalog, but we must then 
apply additional criteria to assess the reliability of each 
detection given that such a catalog will contain many spu- 
rious sources. More conservatively, we choose the low- 
est possible threshold for which contamination by noise is 
unimportant. We aim to produce a catalog with reliable 
colors suitable for robustly modelling of the intrinsic spec- 
tral energy distribution. Using SExtractor version 2.2.2 
(Bertin & Arnouts 1996) with a detection procedure that 
optimizes sensitivity for point-like sources, we construct 
a ii's-band selected catalog with seven band optical-to- 
infrared photometry. 

5.1. Detection 

To detect objects with SExtractor using a constant 
signal-to-noise criterion over the entire image, including 
the shallower outer parts, we divide the point source opti- 
mized /Cg-image by the square root of the weight (exposure 
time) map to create a noise-equalized detection image. A 
source enters the catalog if, after low-pass filtering of the 
detection image, at least one pixel is above « 5 times the 
standard deviation of the filtered background, correspond- 
ing to a total /^s-band magnitude limit for point sources 
of Kg K, 26.0. This depth is reached for the central 4.5 
arcmin^. In total we have 833 detections in the entire sur- 
vey area of 8.3 arcmin^. Initially 820 sources are found 
but the detection software fails to detect sources lying in 
the extended wings of the brightest objects. To include 
these, we fit the surface brightness profiles of the brightest 
sources with the GALPHOT package (Pranx et al. 1989) 
in IRAF, subtract the fit, and carry out a second detection 
pass with identical parameters. Thirteen new objects enter 
the catalog, and 9 sources detected in the first pass are re- 
placed with improved photometry. The catalog identifica- 
tion numbers of all second-pass objects start at 10001, and 
the original entries of the updated sources are removed. 

Filtering affects only the detection process and the 
isophotal parameters; other output parameters are af- 
fected only indirectly through bary center and object ex- 
tent. We chose a simple two-dimensional gaussian de- 
tection filter (FWHM= 0'.'46), approximating the core of 
the effective Xg-band PSF well. Hence, we optimize de- 
tectability for point-like sources, introducing a small bias 
against faint extended objects. In principle it is possible 
to combine multiple catalogs created with different filter 



sizes but merging these catalogs consistently is a compli- 
cated and subjective process yielding a modest gain only 
in sensitivity for larger objects. We prefer the small fil- 
ter size equal to the PSF in the detection map because 
the majority of faint sources that we detect are compact 
or unresolved in the NIR and because we wish to mini- 
mize the blending effect of filtering on the isophotal pa- 
rameters and on the confusion of sources. SExtractor ap- 
plies a multi-thresholding technique to separate overlap- 
ping sources based on the distribution of the filtered Ks- 
band light. About 20% of the sources are blended because 
of the low value of the isophotal threshold; in SExtractor 
this must always equal the detection threshold. With the 
deblending parameters we used, the algorithm succeeds in 
splitting close groups of separate galaxies, without "over- 
splitting" galaxies with rich internal structure. 

We tested sensitivity to false detections by running SEx- 
tractor on a specially constructed Xs-band noise map cre- 
ated by subtracting, in pairs, individual Xg-band images 
of comparable seeing, after zero point scaling, and coav- 
eraging the weighted difference images. This noise im- 
age has properties very similar to the noise in the original 
reduced image, including contributions from the detector 
and reduction process, but with no trace of astronomical 
sources. Our detection algorithm resulted in only 11 spu- 
rious sources over the full area. 



5.2. Optical and NIR Photometry 

We use SExtractor's dual-image mode for spatially accu- 
rate and consistent photometry, where objects are detected 
and isophotal parameters are determined from the K^- 
band detection image while the fluxes in all seven bands 
are measured in the registered and PSF matched images. 
We used fluxes measured in circular apertures APER{D) 
with flxed-diameters D, isophotal apertures APER{ISO) 
determined by the i^g-band detection isophote at the 
5(7 detection threshold, and APER{AUTO) (autoscaling) 
apertures inspired by Kron (1980), which scales an ellip- 
tical aperture based on the first moments of the Ks-hand 
light distribution. We select for each object the best aper- 
ture based on simple criteria to enable detailed control of 
photometry. We define two types of measurements: 

• "color" fiux, to obtain consistent and accurate 
colors. The optimal aperture is chosen based on 
the Kg flux distribution, and this aperture is used 
to measure the flux in all other bands. 

• "total" flux, only in the ifg-band, which gives the 
best estimate of the total Kg flux. 



For both measurements wc treat blended sources differ- 
ently from unblended ones, and consider a source blended 
when its BLENDED or BIAS flag is set by SExtractor, as 
described in Bertin & Arnouts (1996). 

Our color aperture is chosen as follows, introducing 
the equivalent of a circular isophotal diameter Di^o = 
'2'\J Aisoji^ based on Aiso, the measured non-circular 
isophotal area within the detection-isophote: 
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if the source is not blended 

(APER{ISO) (0'/7 < Diso < 2'.'0) 
APER{COLOR) = I APER (0'/7) (Aso < 0'/7) 
[aPER (2'/0) (Aso > 2'/0) 

else if the source is blended 

(APER (Diso/s) {(f.'l < Di,o/s < 
APER{COLOR) = i APER (0'/7) {Di,o/s < 0'/7) 
[aPER (2'/0) (A«o/s > 2'.'0) 

(4) 

The parameter s is the factor with which we shrink the 
circular apertures centered on blended sources, increasing 
the separation to the blended neighbour such that mutual 
flux contamination is minimal. This factor depends on 
the data set, and for our ISAAC Kg image we find that 
s = 1.4 is most successful. The smallest aperture con- 
sidered, APER(0'.'7), w 1.5 FWHM of the effective PSF, 
optimizes the S/N for photometry of point sources in un- 
weighted apertures and prevents smaller more error-prone 
apertures. The largest allowed aperture, APER(2'.'0), pre- 
vents large and inaccurate isophotal apertures driven by 
the filtered Ks-light distribution. We continuously as- 
sessed the robustness and quality of color flux measure- 
ments by inspecting the fits of redshifted galaxy templates 
to the flux points, as described in detail in §6. 

We calculate the total flux in the Ks-hand from the flux 
measured in the AUTO aperture. We define a circularized 
AUTO diameter Dauto = 2-\/^auto/7r with Aauto the area 
of the AUTO aperture, and define the total magnitude as: 

if the source is not blended 



To physically interpret the seven-band photometry for 
our Ks-ha.nd selected sample, we use a photometric red- 
shift {Zphot) technique explained in detail by ROl. In sum- 
mary, we correct the observed flux points for galactic ex- 
tinction (see Schlegel, Finkbeiner, & Davis 1998) and we 
model the rest-frame colors of the galaxies by fitting a lin- 
ear combination of redshifted empirical galaxy templates. 
The redshift with the lowest statistic, where 



2'.'0) 
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(6) 



APER{TOTAL) = APER{AUTO) 
else if the source is blended 
APER{TOTAL) = APER{COLOR) 



(5) 



Finally, we apply an aperture correction using the 
growth curve of brighter stars to correct for the flux lost 
because it fell outside the "total" aperture. This aper- 
ture correction is necessary because it is substantial for 
our faintest sources, as shown in Figure 5 where we com- 
pare different methods to estimate magnitude. The aper- 
ture correction reaches 0.7 mag at the faint end, therefore 
magnitudes are seriously underestimated if the aperture 
correction is ignored. 

We derive the la photometric error for all measurements 
from Eq. 3 with the best-fit values shown in Table 3. 
These errors may overestimate the uncertainty in colors 
of adjacent bands (see section 4.4) but it should represent 
well the photometric error over the entire 0.3/i — 2.2y^m 
wavelength range. The magnitudes may suffer from ad- 
ditional uncertainties because of surface brightness biases 
or possible biases in the sky subtraction procedure which 
could depend on object magnitude and size. 

6. PHOTOMETRIC REDSHIFTS 



is then chosen as the most likely Zphot- Using a linear 
combination of SEDs as F'^"'^"^ minimizes the a priori as- 
sumptions about the nature and stellar composition of the 
detected sources. 

Our data set with three deep NIR bands samples the 
position of Balmer/4000A break over 1 ^ z < 4, allow- 
ing us to probe the redshift distribution of more evolved 
galaxy types that may have little rest-frame UV flux and 
hence a weak or virtually absent Lyman break. 

6.1. Photometric Templates 

We used the local Hubble type templates E, Sbc, Scd, 
and Im from Coleman, Wu, & Weedman (1980), the two 
starburst templates, SBl and SB2, with low derived red- 
dening from Kinney et al. (1996), and a 10 Myr old sin- 
gle age template model from Bruzual & Chariot (2002). 
The starburst templates are needed because many galax- 
ies even in the nearby Universe have bluer colors than 
the bluest CWW templates. The observed templates are 
extended beyond their published wavelengths into the far- 
ultraviolet by power law extrapolation and into the NIR 
using stellar population synthesis models from Bruzual & 
Chariot (2002), with the initial mass functions and star- 
formation timescales for each template Hubble type from 
Pozzetti, Bruzual, & Zamorani (1996). We accounted for 
internal hydrogen absorption of each galaxy by setting the 
flux blueward of the 912A Lyman limit to zero, and for 
the redshift-dependent cosmic mean opacity due to neu- 
tral intergalactic hydrogen by following the prescriptions 
of Madau (1995). 

6.2. Zphot Uncertainties 

The best test of photometric redshifts is direct compar- 
ison to spectroscopic redshifts, but spectroscopic redshifts 
in the HDF-S are still scarce. We calculate the uncer- 
tainty in the photometric redshift due to the flux mea- 
surement errors using a Monte-Carlo (MC) technique de- 
rived from that used in ROl and fully explained in Rud- 
nick et al. (2002b). At bright magnitudes template mis- 
match dominates the errors, something that is not mod- 
eled by the MC simulation. Hence, the MC error bars for 
bright galaxies are severe underestimates. At fainter mag- 
nitudes, the uncertainty is driven by errors in photome- 
try (Fernandez-Soto, Lanzetta, & Yahil 1999) and the MC 
technique should provide accurate Zphot uncertainties. Ex- 
perience from ROl showed that two ways to correct for the 
template mismatch, setting a minimum fractional flux er- 
ror or setting a minimum Zphot error based on the mean 
disagreement with Zgpec-, either degrade the accuracy of the 
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Fig. 5. — Comparison of methods to estimate total Ks-hand magnitude. Shown are isophotal (top), SExtractor's auto-scaUng AUTO 
(middle), and our "total" magnitudes {bottom) as defined in Eq. 5 and which are aperture corrected using stellar growth curve analysis. We 
subtracted the aperture corrected magnitude measured in an aperture of Cf.'T {M AGAPE R(0'!T) — 0.7), which produces the correct total 
magnitudes for stars and pointlike sources. Stars are marked by star symbols and fluxes are plotted with ±1ct error bars. The turn-up at 
Ks Ki 24 of isophotal and at Kg 25 of the AUTO magnitudes shows that these photometric schemes systematically underestimate the total 
flux at faint levels, due to the decreasing size of the used aperture with magnitude. This effect is nearly absent in the bottom panel, which 
shows the total magnitudes measured in this paper. 



Zphot measurement or reflect the systematic error only in 
the mean, while template mismatch can be a strong func- 
tion of SED shape and redshift. A method based com- 
pletely on Monte-Carlo techniques is preferable because 
it has a straightforwardly computable redshift probabil- 
ity function. This approach is desirable for estimating the 
rest-frame luminosities and colors (Rudnick et al. 2002b). 

Therefore, we modify the MC errors directly using the 
FIRES photometry. In summary, wc estimate the system- 
atic component of the Zphot uncertainty by scaling up all 
the photometric errors for a given galaxy with a constant 
to bring residuals of the fit in agreement with the errors. 
This will not change the best fit redshift and SED and will 
not modify the MC errorbars of faint objects, but it will 
enlarge the redshift interval over which the templates can 
satisfactorily fit the bright objects. Only in case of widely 
different photometric errors between the visible and in- 
frared might the modified MC uncertainties still underes- 
timate the true ZpUot uncertainty. 

In Figure 6 we show a direct test of photometric redshift s 
of the 39 objects in the HDF-S with available spectroscopy 
and good photometry in all bands. The current set of spec- 
troscopic redshifts in the HDF-S will appear in Rudnick et 
al. (2002a). For the small sample that we can directly com- 



pare, we find excellent agreement with no failures and with 
a mean Az/(1 + Zspec) « 0.08 with Az = \zspec - Zphot\- 
It is encouraging to see that the modified 68% error bars 
that were derived from the Monte Carlo simulations are 
consistent with the measured disagreement between Zphot 
and the Zgpec in the HDF-S. However, large asymmetric 
uncertainties remain for some objects, clearly showing the 
presence of a second photometric redshift solution of com- 
parable likelihood at a vastly different redshift, revealing 
limits on the applicability of the photometric redshift tech- 
nique. 

6.3. Stars 

In a pencil beam survey at high Galactic latitude such 
as the HDF-S, a limited number of foreground stars are 
expected. We identify stars as those objects which have 
a better raw foi' single stellar template fit than 
the for the galaxy template combination. The stel- 
lar templates are the NEXTGEN model atmospheres from 
Hauschildt et al. (1999) for main sequence stars with tem- 
peratures of 3000 to lOOOOK, assuming local thermody- 
namical equilibrium (LTE). Models of cooler and hotter 
stars cannot be included because non-LTE effects are im- 
portant. We checked the resulting list of stars using the 
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Fig. 6. — Direct comparison of photometric rcdshifts to the 39 spectroscopic redshifts of objects in the HDF-S with good photometry in 
all bands. The 68% error bars are derived from our Monte Carlo simulations and the diagonal line corresponds to a one-to-one relation to 
guide the eye. While the agreement is excellent with no failures for this small sample and with mean Az/(1 + z) = 0.08, large asymmetric 
uncertainties remain for some objects indicating the presence of a second photometric redshift solution of comparable likelihood at a different 
redshift. 



FWHM in the original i?45o-band image and the Jg — Kg 
color, excluding two objects (catalog IDs 207 and 296) that 
were obviously extended in i?45o, and wc find a total of 57 
stars. As shown in Figure 7, most galaxies are clearly sep- 
arated from the stellar locus in I^u — Kg versus Jg — Kg 
color-color space. Other cooler stars might still resem- 
ble SEDs of redshifted compact galaxies but the latter are 
generally redder in the infrared Jg — Kg than most known 
M or methane dwarfs. Known cool L-dwarfs fall along a 
redder extension of the track traced by M-dwarfs in color- 
color space and have progressively redder Jg — Kg colors 
for later spectral types. However measurements by Kirk- 
patrick et al. (2000) show the L-dwarf sequence abruptly 
stopping a.t{Jg — Kg)jKi2.1 (the subscript noting John- 
son magnitudes, see section 3.1 for the transformations to 
the AB system) whereas even cooler T-dwarfs have much 
bluer {Jg — Kg) J « colors than expected from their tem- 
peratures due to strong molecular absorption. This is im- 
portant because if we would apply a,{Js—Kg)j>2.3 pho- 
tometric criterion to select z > 2 galaxies (as discussed in 
section 8.2), then we should ensure that cool Galactic stars 
are not expected in such a sample. The published data 
on the lowest-mass stars suggest that they arc too blue 
in infrared colors to be selected this way. Only heavily 
reddened stars with thick circumstellar dust shells, such 
as extreme carbon stars or Mira variables, or extremely 
metal-free stars having a hypothetical < 1500K blackbody 
spectrum could also have red {Jg — Kg)j > 2.3 colors but 

http: / / www. strw. leidenuni V. nl / "fires 



it seems unlikely that the tiny field of the HDF-S would 
contain such unusual sources. 

7. CATALOG PARAMETERS 

The ii's-selected catalog of sources is published electron- 
ically. We describe here a subset of the photometry con- 
taining the most important parameters. The catalog with 
full photometry and explanation can be obtained from the 
FIRES homepage 

• ID. — A running identification number in catalog 
order as reported by SExtractor. Sources added 
in the second detection pass have numbers higher 
than 10000. 

• x,y. The pixel positions of the objects 
corresponding to the coordinate system of the 
original (unblocked) WFPC2 version 2 images. 

• RA, DEC. - The right ascension and declination 
in equinox J2000.0 coordinates of which only 
the minutes and seconds of right ascension, and 
negative arcminutes and arcseconds of declination 
are given. To these must be added 22'' (R. A.) and 
-60° (DEC). 

• fcoi,i i cTi- — The sum of counts 

in the "color" aperture fcoi.i in band 

i = {U300, Bi5o,V606, hii, Js, H, Kg} and its 



14 



Labbe et al. 



I 



+ 



1 



(l814- K, 



Fic;. 7. Js — Ks versus Isi4 — Kg color-color diagram (on the AB system) for the sources with Ks < 24 in the HDF-S with a minimum of 
20% of the total exposure time in all bands. Identified stars are marked by a star symbol. The colors are plotted with ±1<t error bars. There 
is a large variation in both / — Ks and Js — Kg colors. Redshifted galaxies are well separated from the stellar locus in color-color space. 



simulated uncertainty Ui, as described in § 5.2. 
The fluxes are given in units of 10~^^ ergs s~^ 
Hz-i cm-2. 

• Ktot ± o-{Ktot)- — Estimate of the total JCg-band 
flux and its uncertainty. The sum of counts in 
the "total" aperture is corrected for missing flux 
assuming a PSF profile outside the aperture, as 
described in § 5.2. 

• ap_col. — An integer encoding the aperture type 
that was used to measure fcoi,i- This is either 

a (1) 0'.'7 diameter circular aperture, (2) 2'.'0 
diameter circular aperture, (3) isophotal aperture 
determined by the detection-image isophote, or 
a (4) circular aperture with a reduced isophotal 
diameter D = sj (yljso/7r)/1.4. 

• apJot. — An integer encoding the aperture type 

that was used to measure Ktot- This is either a 
(1) automatic Kron-like aperture, or a (2) circular 
aperture within a reduced isophotal diameter. 

• fcoUftot- — Circularized radii r = ^J A/tt, 
corresponding to the area A of the specified "color" 
or "total" aperture. 

• Aigo, Aauto- — Area of the detection isophote 
Aiso and area of the autoscaling elliptical aperture 
Aauto = TT * a* b with semi-major axis a and 
semi-minor axis b. 



• FWHMk,FWHMi. — Full width at half 
maximum of a source in the Ks detection image 
FWHMk, and that of the brightest /814-band 
source that lies in its detection isophote FWRMj. 
We obtained the latter by running SExtractor 
separately on the original J8i4-image and cross 
correlating the /814-selected catalog with the 
ifg-limited catalog. 

• Wi. — The weight Wi represents, for each band 
i, the fraction of the total exposure time at the 
location of a source. 

• flags.-- Three binary flags arc given. The bias 
flag indicates either that the AUTO aperture 
measument is affected by nearby sources, or marks 
apertures containing more than 10% bad pixels. 
The blended flag indicates overlapping sources, 
while the star flag shows that the source SED is 
best fit with a stellar template (see section 6.3). 

8. ANALYSIS 

8.1. Completeness and Number Counts 

The completeness curves for point-sources in the Js and 
Ks-hsLud as a function of input magnitude are shown in 
Figure 8. Our 90% and 50% completeness levels on the 
AB magnitude system are 25.65 and 26.25, respectively, 
in Ks, and 26.30 and 26.90 in J^. 
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Fig. 8. — Completeness curves (on the AB system) for the detectability of point sources in Js (triangles) and Kg (points), based on 
simulations where we calculated the recovered fraction of stars that were dimmed to magnitudes between 22 and 28 and embedded in the 
survey images. The detection threshold of the source extraction software was set to 3.5a of the filtered background rms. The dotted lines 
indicate the 50% and 90% completeness levels. 



We derived the limits from simulations where we ex- 
tracted a bright non-saturated star from the survey image 
and add it back 30000 times at random locations, apply- 
ing a random flux scaling drawn from a rising count slope 
(or an increasing surface density of galaxies with magni- 
tude) to bring it to magnitudes between 22 < Kg^AB < 28. 
We added the dimmed stars back in series of 30 realiza- 
tions so that they do not overlap each other. The rising 
count slope needs to be considered because the slope in- 
fluences the number of recovered galaxies per apparent 
magnitude, as described below. The input count slope is 
based on the observed surface densities in the faint mag- 
nitude range where the signal-to-noise is 60 < SNR < 10 
(or 23 < Ks^tot ^ 25) and where incompleteness docs not 
yet play a role. We used only the deepest central 4.5 
arcmin^ of the Jg and Ks images {w > 0.95) with near uni- 
form image quality and exclude four small regions around 
the brightest stars. In the simulation images we extract 
sources following the same procedures as described in §5.1, 
but applying a reduced (« 3.5cr) detection threshold. We 
measure the recovered fraction of input sources against ap- 
parent magnitude, and from this we estimate the detection 
efficiency of point-like sources which we use to correct the 
observed number counts. We executed this procedure in 
the Js and i^g-band. 

The resulting completeness curves assume that the true 
profile of the source is point-like and therefore they should 
be considered upper limits. An extended source would 
have brighter completeness limits depending on the true 
source size, its flux profile, and the filter that is used in the 
detection process. However, the detailed treatment of de- 
tection efiiciency as a function of source morphology and 



detection criteria is beyond the scope of this paper. 

When using the completeness simulations to correct the 
number counts, we choose a simple approach and apply a 
single correction down to « 50% completeness, based on 
the ratio of the simulated counts per input magnitude bin 
to the total recovered counts per observed total magnitude 
bin. More sophisticated modeling is possible but requires 
detailed knowledge of the intrinsic size and shape distribu- 
tion of faint NIR galaxies. The simple approach corrects 
for all effects resulting from detection criteria, photomet- 
ric scheme, incompleteness, and noise peaks. We find it 
works well if the total magnitude of sources is measured 
correctly, with little systematic difference between the in- 
put and recovered magnitudes, which is the case for our 
photometric scheme (see section 5.2). It is worth noticing 
that at Ks,AB > 25 we actually recovered slightly more 
counts in the observed magnitude bins than we put in. 
This is caused by the fact that, in the case of a rising 
count slope, there are more faint galaxies boosted by posi- 
tive noise peaks than bright galaxies lost on negative noise 
peaks. This effect is strong at low signal-to-noise fluxes 
and results in a slight excess of recovered counts. This is 
the main reason that we required little correction up to 
the detection threshold, except for the faintest 0.5 mag 
bin centered on Kg = 26.0, which contained false positive 
detections due to noise. After removing stars (see section 
6.3), we plot in Figure 9 the raw and corrected source 
counts against total magnitude. 

Figure 10 presents a compilation of other deep K-hand 
number counts from a number of published studies. The 
FIRES counts follow a d/ogf(N)/dm relation with a loga- 
rithmic slope a w 0.25 at 20 < Ks,j < 22 (Johnson mag- 
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Fig. 9. — Differential Ks-ha,nd counts (on tlic AB system) of galaxies in the HDF-S. The counts are based on auto-scaling apertures (Kron 
1980) for isolated sources and adapted isophotal apertures i^or blended sources, both corrected to total magnitudes using stellar growth curve 
measurements. Raw counts {open circles) and counts corrected for incompleteness and false positive detections using point source simulations 
{filled circles) are shown. The small corrections at magnitudes > 23 reflect missed sources due to confusion. Effective corrections at the 
faintest magnitudes Ks ^ 25 — 25.5 are very small because the loss of sources on negative noise regions (incompletenesss) is compensated by 
the number of sources pushed above the detection limit by positive noise fluctuations. Only the faintest 0.5 mag bin centered on Ks = 26.0, 
bordering the 3.5a detection limit {Kg « 26.3), is significantly corrected because of a contribution of false positive detections. 



nitudes) and decline at fainter magnitudes to a « 0.15 at 
22.0 < Ks,j < 24. This flattening of the slope has not 
been seen in other deep NIR surveys, where we emphasize 
that the FIRES HDFS field is the largest and the deepest 
amongst these surveys, and that only the counts in the last 
FIRES bin at Ks,j = 24.25 were substantially corrected. 
It is remarkable that the SUBARU Deep Field count slope 
a ~ 0.23 of Maihara et al. (2001) looks smooth compared 
to the HDF-S although their survey area and the raw count 
statistics are slightly smaller. 

Other authors (Djorgovski et al. 1995; Moustakas et al. 
1997; Bershady, Lowenthal, & Koo 1998) find logarith- 
mic counts slopes in K ranging from 0.23 to 0.36 over 
20 < Kj < 23 — 24, however the counts in the faintest bins 
in these surveys were boosted by factors of 2 — 3, based 
on completeness simulations. The origin of the faint-end 
discrepancies of the K counts is unclear. Cosmic vari- 
ance can play a role, because the survey areas never ex- 
ceed a few armin^, but also differences in the used fil- 
ters {Kg, Kf, K) and differences in the techniques and as- 
sumptions used to estimate the total magnitude (see 5) 
or to correct the counts for incompleteness may be im- 
portant. Further analysis is needed to ascertain whether 
size-dependent biases in the completeness correction play 
a role in the faint-end count slope. 

8.2. Color-Magnitude Distributions 



Figures 11-14 show color-magnitude diagrams of Kg- 
selected galaxies in the HDF-S. The /814 — Kg versus Kg 
color-magnitude diagram in Figure 11 shows a large num- 
ber of extremely red objects (EROs) with Ign — Kg > 2.6 
(on the AB system) or (/814 — Ks)j > 4 (Johnson). 
There appears to be an excess of EROs at total magni- 
tudes Kg^AB ~ 23 compared to magnitudes Kg_AB ~ 24. 
This is not caused by the insufficient signal-to-noise ra- 
tio in the /814 measurements. In a similar diagram for 
the Jg — Kg color shown in Figure 12, there is a strik- 
ing presence at the same Kg magnitudes of sources with 
very red {Jg — Kg)AB ^ 1-34 or [Jg — Kg)j > 2.3 colors. 
Such sources were also found by Saracco et al. (2001), us- 
ing shallow NIR data, who suggested they might be dusty 
starbursts or ellipticals at z > 2. Interestingly, any evolved 
galaxy with a prominent Balmer/4000A discontinuity in 
their spectrum, like most present-day Hubble Type galax- 
ies, would have such very red observed NIR colors if placed 
at redshifts z > 2. While the {Jg — Kg)j > 2.3 sources we 
find are generally morphologically compact, with excep- 
tions, we do not expect the sources with good photometry 
to be faint cool L-dwarf stars because known colors of such 
stars are (J, — Kg),} < 2.1 (see section 6.3). The photo- 
metric redshifts of all red NIR galaxies are Zphot ^ 2, but 
they would be missed by ultraviolet-optical color selection 
techniques such as the U-dropout method, because most of 
them are barely detectable even in the deepest optical im- 
ages. One bright NIR galaxy is completely undetected in 
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Fig. 10. — FIRES Xs— band galaxy counts (on the conventional Johnson system) compared to published counts in deep ii'-band fields. The 
corrected counts {filled circles) are shown for FIRES data. The Maihara ct al. (2001) counts have been plotted to their S/N^ 3 limit. The slope 
at magnitudes Kj > 21 is flatter than reported in other surveys although straightforward comparisons are difficult, due to model-dependent 
correction factors of ^ 2 — 3 applied to the faintest data points in these surveys. The nature of the scatter in count slopes is unclear but 
field-to-field variations as well as different photometry and corrections procedures likely play a role. The FIRES counts need little correction 
for completeness effects or false positive detections, except for the K^^j = 24.25 bin. 



the original WFPC2 images. The {Js — Ks)j > 2.3 sources 
are studied in more detail by Franx et al. (2002) and the 
relative contributions of these galaxies and U-dropouts to 
the rest-frame optical luminosity density will be presented 
in Rudnick et al. (2002b). 

If we select sources with 1.95 < Zphot < 3.5, we find 
clear differences in the Veoe — H versus H color-magnitude 
diagram between our NIR-selected galaxies in the HDF-S 
and those of the HDF-N (compare Figure 14 to Figure 1 
of Papovich, Dickinson, & Ferguson 2001). Over a simi- 
lar survey area and to similar limiting depths, we find 7 
galaxies redder than (Vgoe — H)ab ^3 and brighter than 
total magnitude Hab S 25.5, compared to only one in the 
HDF-N. 

While the surface density of such galaxies is not well 
known, it is clear that the HDF-N contains far fewer of 
them than the HDF-S. It remains to be seen if this is 
just field-to-field variation, or that one of the two fields 
is atypical. The results of the second much larger FIRES 
field centered on MS1054-03 (Forster Schreiber et al. 2002) 
should provide more insight into this issue. We note that 
all 7 (Veoe ^ H)ab ^ 3 galaxies in the HDF-S are also 
amongst the brightest 16 [Jg — Ks).j > 2.3 sources. Fig- 
ure 15 shows the surface densities of EROs and galaxies 
with [Js — Kg) J > 2.3 colors as function of iiTg-band total 
magnitude. The surface density of EROs peaks around 
Kg.AB ~ 23 and then drops or flattens at fainter magni- 
tudes, contrary to the number of{Js — Ks)j > 2.3 galaxies 
which keeps rising to the faintest Kg total magnitudes in 



our catalog. 

9. SUMMARY AND CONCLUSIONS 

We have presented the results of the FIRES deep NIR 

imaging of the WFPC2-field of the HDF-S obtained with 
ISAAC at the VLT: the deepest ground-based NIR data 
available, and the deepest Ks-hamd of any field. We con- 
structed a i^s-sclected mtilticolor catalog of galaxies, con- 
sisting of 833 objects with Ks^ab ^ 26 and photometry 
in seven-bands from 0.3 to 2.2/z for 624 of them. These 
data are available electronically together with photomet- 
ric redshifts for 567 galaxies. Our unique combination of 
deep optical space-based data from the HST together with 
deep ground-based NIR data from the VLT allows us to 
sample light redder than the rest-frame V-band in galaxies 
with z < 3 and to select galaxies from their rest-frame op- 
tical properties, obtaining a more complete census of the 
stellar mass in the high-redshift universe. We summarize 
our main findings below; 

• The Xg-band galaxy counts in HDF-S turn over at 
the faintest magnitudes and fiatten from a « 0.25 
at total AB magnitudes 22 < Jf^ < 24 to a w 0.15 
at 24 < JCs < 26; this is flatter than counts in 
previously published deep NIR surveys, where the 
FIRES HDF-S field is largest and deepest amongst 
these surveys. The nature of the scatter in the 
faint-end counts is yet unclear but field-to-field 
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Fig. 11. — 7814 — versus Kg color-magnitude relation (on the AB system) for Xs-selected objects in the HDF-S. Only sources with a 
minimum of 20% of the total exposure time in all bands are included and identified stars are marked by a star symbol. Colors are plotted with 
±1(7 error bars, and Isi4 measurements with S/N < 2 (triangles) are plotted at their 2a confidence interval, indicating lower limits for the 
colors. There are more red sources with /gi4 — Kg > 2.6 at ~ 23 than at at ii" ^ 24 where the /814 is still sufficiently deep to select them. 
The transformation of the 7814 — color from the AB system to the Johnson magnitude system is (78i4 — Ks)j = {I814 — Ks)ab + 1-43. 
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Fig. 12. — Same as Figure 11 for the Js — Ks color. Striking is the number the galaxies with very red NIR colors Js — Kg > 1.34 (on 
the AB system) or Jg — Ks > 2.3 (Johnson). These systems have photometric redshifts z > 2 and are extremely faint in the observer's 
optical; as such they would not be selected with the U-dropout technique. Identified stars are well separated from redshifted galaxies and 
almost all have Js — Ks < colors. The transformation of the Js — Ks color from the AB system to the Johnson magnitude system is 
[Js - Ks)j = {Js - Ks)ab + 0.96. 
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Fig. 13. — Same as Figure 11 for the H — Kg color. One of the galaxies is extremely red with H — Kg 2.2 and is barely visible in Js and 
H. The transformation of the H — Kg color from the AB system to the Johnson magnitude system is {H — Ks)j = {H — Ks)ab + 0.48. 
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Fig. 14. — Veo6 — H versus H color-magnitude diagram (on the AB system) for galaxies in the HDF-S Xs-selected catalog with 
1.95 < ^p^ot < 3.5. Filled symbols indicate galaxies with spectroscopy. The number of candidates for red, evolved galaxies is much 
higher than in the HDF-N for a similar survey area, as shown in a identical plot in Fig. 1 of Papovich, Dickinson, & Ferguson (2001): we 
find 7 galaxies redder than Vgoe AB ~ > 3 and brighter than Hab ^ 25.5, compared to only one in the HDF-N. Galaxies with S/N 

< 2 for the Veoe measurement (triangles) are plotted at the 2<t confidence limit in Vgog, indicating a lower limit on the Veoe — H color. The 
subsample of galaxies having red (Js — Ks)j > 2.3 colors {open squares) is also shown. The transformation of the Veoe — Hs color from the 
AB system to the Johnson magnitude system is (Veoe — H)j = (Veoe — H)ab + 1.26. 



22 



Labbe et al. 



10.0 



1.0 



0.1 



18 20 22 24 26 

Kg total magnitude 



Fig. 15. — The surface densities of galaxies selected by color in the 5cr catalog of sources of HDF-S. Presented are galaxies with 
{Js — Ks)j > 2.3 (in Johnson magnitudes) {diamonds), extremely red objects with (/814 — Ks)j > 4 (squares), and all Xs-selccted galaxies 
{filled points) as a function total /sTs-band AB magnitude. Only sources with a minimum exposure time of 40%, 40% and 90% of the total 
in /814, Js and Ks arc plotted, so that the selection in Kg is uniform over the area, and the /gi4 and Js observations arc sufficiently deep to 
prevent a bias against objects with very red 7814 — Ks and Js — Ks colors. No corrections to the counts have been applied. The errorbars are 
poissonian and might underestimate the true uncertainty which would also contain contributions from large scale structure. 



varations as well as different analysis techniques 
likely play a role. 

• The HDF-S contains 7 sources redder than 
(^06 — H)ab ^ 3 and brighter than total 
magnitude Hab ^ 25.5 at photometric redshifts 
1.95 < Zphot < 3.5, while such galaxies were 
virtually absent in the HDF-N. They are much 
redder than regular U-dropout galaxies in the same 
field and are candidates for relatively massive, 
evolved systems at high redshift. The difference 
with the HDF-N might just reflect field- to-field 
variance, calling for more observations to similar 
limits with full optical-to-infrared coverage. 
Results from the second and larger FIRES field 
centered on MS1054-03 (Forster Schreibcr et al. 
2002) should provide more insight into this issue. 

• We find substantial numbers of red galaxies with 
{Js — Kg) J > 2.3 that have photometric redshifts 
Zphot > 2. These galaxies would be missed by 
ultraviolet-optical color selection techniques such 
as the U-dropout method because most of them 
are barely detectable even in the deepest optical 
images. The surface densities of these sotirces in 
our field keeps rising down to the detection limit 
in Kg, in contrast to the number counts of EROs 



which peak at Ks.ab ~ 23 and then drop or flatten 
at fainter magnitudes. 

The results of the HDF-S presented in this paper demon- 
strate the necessity of extending optical observations to 
near-IR wavelengths for a more complete census of the 
early universe. Our deep i^s-band data prove invaluable 
for they probe well into the rest-frame optical at 2 < z < 4, 
where long-lived stars may dominate the light of galaxies. 
We are pursuing follow-up programs to obtain more spec- 
troscopic redshifts needc;d to confirm the above results. 
Updates on the FIRES programme and access to the re- 
duced images and catalogues can be found at our website 
http : / / www. st r w . leidenuniv . nl/ ~ fir es . 

The data here presented have been obtained as part of 
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motis efforts in taking these data and making them avail- 
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their dedicated work. This work was supported by a grant 
from the Netherlands Organization for Scientific Research. 
We would like to thank the Lorentz Center of Leiden Uni- 
versity for its hospitality during several workshops. The 
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Labbe et al. 



Table 1 



Summary of the HDF-S Observations 



Camera 


Filter 


Number of 


Integration Time 


Image Quality^ 






Exposures 


(h) 


(arcsec) 


WFPC2 


F300W 


102 


36.8 


o'.'ie 


WFPC2 


FA50W 


51 


28.3 


0'.'14 


WFPC2 


F606W 


49 


27.0 


0'.'13 


WFPC2 


FSIAW 


56 


31.2 


0'.'14 


ISAAC 


Js 


1007 


33.6 


0'.'45 


ISAAC 


H 


968 


32.3 


0'.'48 


ISAAC 


Ks 


2136 


35.6 


0'.'46 



The full width at half-maximum of the best-fitting Gaussian. 



Table 2 
Zero Points for the HDF-S 



Zero Point 
Data Set Johnson mag AB mag 





19.43 


20.77 


-B450 


22.01 


21.93 




22.90 


23.02 




21.66 


22.09 


Js 


24.70 


25.60 


H 


24.60 


25.98 


Ks 


24.12 


25.98 



Table 3 





BACKCROrKD 


KOISE 


L\ THE 


IIDF-S imaces 


Data Set^ 


rms background variation^ 


a3 


b3 


la sky noise limit* 


C^300 


1.34e-05 


2.51 


0.38 


29.5 


-B450 


1.88e-05 


2.65 


0.40 


30.3 




3.80e-05 


2.49 


0.39 


30.6 




2.79e-05 


2.59 


0.39 


30.0 


Js 


0.0069 


1.46 


0.047 


28.6 


H 


0.0165 


1.43 


0.044 


28.1 


Ks 


0.0163 


1.49 


0.038 


28.1 



^The images are all at the 0'.'119 pixel ^ scale. The WFPC2 are 3x3 block summed, and all 
images are smoothed to match the image quality of the //-band. 

^Pixcl-to-pixcl rms variations (in instrumental counts per second) as measured directly in empty 
parts of the registered convolved images which were used for photometry. 

^Best-fit parameters of Eq. 3 which gives the effective rms variation of the background as a 
function of linear size of the aperture. 

*The lc7 sky noise limit in a 0'.'7 circular diameter aperture (« 0.4 arcsec^) in AB magnitudes 
using Eq. 3. 



